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We propose a theory of adiabaticity in quantum Markovian dynamics based on a decomposition of 
the Hilbert space induced by the asymptotic behavior of the Lindblad semigroup. A central idea of 
our approach is that the natural generalization of the concept of eigenspace of the Hamiltonian in the 
case of Markovian dynamics is a noiseless subsystem with a minimal noisy cofactor. Unlike previous 
attempts to define adiabaticity for open systems, our approach deals exclusively with physical entities 
and provides a simple, intuitive picture at the underlying Hilbert-space level, linking the notion of 
adiabaticity to the theory of noiseless subsystems. As an application of our theory, we propose a 
framework for decoherence-assisted computation in noiseless codes under general Markovian noise. 
We also formulate a dissipation-driven approach to holonomic computation based on adiabatic 
dragging of subsystems that is generally not achievable by non-dissipative means. 
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Introduction. — The adiabatic theorem is a simple and 
powerful result that has been known since the early days 
of quantum mechanics [H, 0]- It states roughly that a 
closed system in an eigenstate of a continuously per- 
turbed Hamiltonian remains in an instantaneous eigen- 
state in the limit of slow perturbations if the correspond- 
ing eigenvalue is separated from the rest of the spectrum 
by a gap. Adiabaticity in quantum mechanics has ap- 
plications in a wide range of areas, including quantum 
chemistry [3|] , geometric phases 0, Q , quantum Hall ef- 
fect [1], STIRAP 0, and quantum phase transitions Q. 
More recently, the adiabatic theorem has been the sub- 
ject of increased interest in relation to quantum infor- 
mation processing, where it has served as a basis for a 
variety of schemes, including holonomic quantum com- 
putation and adiabatic quantum algorithms [10|]- 

Given the importance of the concept of adiabaticity in 
closed quantum systems, it is natural to ask how this con- 
cept extends to the dynamics of systems interacting with 
an environment. This question is of particular interest 
from the point of view of quantum information processing 
where decoherence is a major obstacle to the construc- 
tion of reliable quantum devices, and at the same time 
non-uni tary processes are an important tool for quantum 
control In Ref. Sarandy and Lidar proposed 

an approach to the adiabatic dynamics of open quantum 
systems, defining adiabaticity as the regime in which the 
operator subspaces corresponding to the instantaneous 
Jordan blocks of the generator of the dynamics evolve 
independently (for adiabaticity in weakly open systems, 
see Ref. [13] )■ This definition is motivated by the formal 
analogy between the Schrodinger equation and the time- 
dependent Markovian master equation written in a coher- 
ence basis, both being first-order linear vector differential 
equations with the difference that the generator of the 
master equation is generally not diagonalizable (hence 
the Jordan decomposition). But while in closed sys- 
tems the phenomenon of adiabaticity concerns the decou- 



pled evolution of eigenspaces of the Hamiltonian which 
themselves are Hilbert spaces containing physical states, 
the Jordan blocks correspond to generally nonorthogonal 
subspaces of the space of linear operators that need not 
contain density matrices or even observables and may de- 
cay to zero even when mutually decoupled. In the present 
paper, we propose a different approach, based primarily 
on physical considerations, which yields an inequivalent 
picture of open-system adiabaticity that links adiabatic 



dynamics to the theory of noiseless subsystems 

Taking as a ground the basic physical characteristic of 
adiabatic closed-system evolutions — namely, that these 
are quasi-static evolutions where under sufficiently slow 
changes of the Hamiltonian a system in a stationary state 
evolves so as to remain in a stationary state with re- 
spect to the changed Hamiltonian — we look for a gener- 
alization of this phenomenon to the case of Markovian 
dynamics. The key insight of our approach is that the 
natural generalization of the eigenspaces of the Hamilto- 
nian corresponding to distinct eigenvalues are noiseless 
subsystems whose noiseful cofactors support unique fixed 
states. A decomposition of the Hilbert space into such 
subsystems arises naturally from the asymptotic behavior 



of the Lindblad semigroup [15[. We define adiabaticity 



as the regime in which the stationary states over such 
a noiseless subsystem and its cofactor remain stationary 
with respect to the Lindbladian as it changes. We derive 
an adiabatic theorem based on this definition. 

To illustrate the utility of our formalism, we propose 
two applications. One is a framework for decoherence- 
assisted computation in noiseless codes which general- 
izes the approach of Beige et al. [31 to subsystems and 
general noise models. The other is a dissipation-driven 
approach to holonomic quantum computation based on 
adiabatic "dragging" of subsystems 17| along paths that 
are generally not achievable by non-dissipative means. 

Generalization of eigenspaces. — Our starting point is 
the observation that the eigenstates of a Hamiltonian H 
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are the stationary state vectors of its dynamics. In par- 
ticular, ah stationary density matrices under the evolu- 
tion dp/dt = —i[H,p\ (we set h = 1) have the direct-sum 
form p = 0jPi/3i, J2iPi = 1) Pi > 0, where pi are den- 
sity matrices over the eigenspaces Hi of H corresponding 
to distinct eigenvalues. In more general quantum pro- 
cesses, the stationary states are organized as operators 
over noiseless subsystems tensored with a fixed de nsity 
matrix over the corresponding noiseful co-subsystem [iSj . 
Consider a time-homogenous finite-dimensional Marko- 
vian dynamics described by the Lindblad equation ,19j 



dp 
Itt 



-J^ = -i[H,p] + 



]{Upl\-\l\up-^-pl\l,)^Cp, 

(1) 



where Li are Lindblad operators. As shown in Ref. [11 
Eq. ([Ij induces a decomposition of the Hilbert space 



(2) 



where T-Lf^ are noiseless subsystems jl4{ , "Hf are noiseful 
subsystems that support unique fixed states, and /C is a 
decaying subspace. More particularly, it was shown that 
for any initial state p(0), the solution of Eq. (HI satisfies 

3{Pk,pt} ■■ lim \p{t) - 0p,-e-^<V/e'<* ® I = 0, 

j 

(3) 

where are density matrices on the unitarily noiseless 
subsystems Tif = ® ^ T-Lf, evolving under the Hamilto- 



of H{s) with multiplicity m, and Pi{s) be the (twice- 
differentiable) projector on the corresponding eigenspace 
Hi{s) — Pi{s)T-L. [Note that m = const(s) implies that 
ti{s) is separated from the rest of the spectrum by a 
gap. The adiabatic theorem has been extended to cases 
without a gap [2l|, but in this paper we restrict to the 
standard formulation.] The eigenspace 'Hi{t/T) is said to 
evolve adiabatically under H{t/T) if any state initially 
in HiiO) remains in Hiit/T), t e [0,r]. Equivalently, if 
we change the basis via a unitary U{s) so that Pi be- 
comes fixed, in the new basis the dynamics is driven by 
the effective Hamiltonian H'{t/T) = H{t/T) + ^V{t/T), 
where H{s) = U{s)H{s)U{s)^ = er{s)P; + H^is) with 
H^{s) having support on the orthogonal complement of 
Hi, and V{s) = i^^^^J7^(s). Adiabaticity then refers to 
the regime in which any state initially in Hi remains in 
Hi despite the action of ^V{t/T). The adiabatic theo- 
rem states f3| that in the limit of large T, one approaches 
perfect adiabaticity where the states in Hi evolve via the 
unitary Ui{s) — Texp(— « PiV{q)Pidq) where T de- 
notes time ordering. The error scales with T as 0( j^), 
where A > is a fixed energy scale (e.g., the energy gap). 

Note that unlike the "folk" adiabatic condition which 
is known to be insufficient [2^ . this theorem (similarly 
to the one derived below) is concerned with the scaling 
of the error as a function of T for a fixed curve H{s). 

Adiabaticity in Markovian dynamics. — Consider a 
time-dependent Lindbladian C{t/T) changing along a dif- 
ferentiable curve £(s), s € [0,1]. For every s, £{s) 
induces a decomposition of the Hilbert space H = 
®ijHfj{s) (g) Hf{s) © /C(s) as explained earlier. Let 
H^iis) and Hf{s) [dimH^i{s) = m, dimHf (s) = n] be 



nians Hj, are fixed full-support states on Hf , and subsystems of the type above, and let Vkiis) \Pki{s)p — 



J2kPk = Pk > 0. The noiseless subsystems Hfj are the 
eigenspaces of Hf. The stationary states have the form 



P=®^^P^3P^ 



^ J^inPtj = > 0, where pf are 



density matrices on Hfj. This suggests that the subsys- 
tems Hfj whose cofactors Hf support unique fixed states 
gf can be thought of as the generalization of eigenspaces 
corresponding to distinct eigenvalues. 

How do we find the decomposition ([2]) for a given Lind- 
bladian C7 An algorithm for finding the noiseless subsys- 
tems of a completely positive trace-preserving (CPTP) 
map that runs in time 0[(dim'H)^] was described in 
Ref. [l^ (see also Ref. [l^]). It is based on finding the 
left and right operator eigenspaces corresponding to the 
eigenvalue 1 of the CPTP map. Since Eq. ^ is equiv- 
alent to the continuous application of an infinitesimal 
CPTP map, the same algorithm can be used here (the 
eigenvalue 1 of the map translates to eigenvalue of £). 

Before we introduce adiabaticity for Markovian dy- 
namics, let us briefiy review the closed-system case. 

Adiabaticity in closed systems. — Consider a time- 
dependent Hamiltonian H{t/T) changing along a differ- 
entiable curve H{s), s € [0, 1]. Let ei{s) be an eigenvalue 



TrB{Pfer(s)pPfc1^(s)}® ef (s) where P,^,^{s) is the pro- 
jector on Hfi{s) ® Hf {s) and Tr^ denotes partial trace 
over Hf] be the (twice-differentiable) superoperator pro- 
jector on the fixed points over Hfi{s) Hf{s). 

Note. Similarly to the closed-system case, the assump- 
tion that dim'H^;(s) and dim'HP(s) are constant implies 
an analogue of the gap condition (see Appendix A). 

Definition. The noiseless subsystem Hfi{t/T) and 
its noisy cofactor Hf{t/T) evolve adiabatically under 
C{t/T), if any state over ^^((0) (g) Hf{Q) of the form 

tate p{t) 

gfit/T) over HUt/T) ®Hf [t/T), t e [0,T]. 

As in the case of closed systems, it is convenient to 
consider a basis rotated by a unitary U{s), in which Hfi 
and Hf are fixed. In this basis, the master equation is 

^=-^[Vit/np]+Cit/T)p, (4) 

where C{s) is the Lindbladian with H(s) replaced 
by U{s)H{s)U{s)'< and L,{s) by U{s)Li{s)U{s)\ and 



p(0) = p{0)fi (g gf{0) evolves to a state p{t) = p(t)^, 



V{s) 



■ dU{s) 



[/(s)t. (We will not use a different no- 



tation for p in this basis but will keep in mind the basis 
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we are working in.) Adiabaticity then means that any 
state p(0) = Pfc/(0) (g) (0) remains of the form p{t) = 
Pkii^) '8> Qf{t/T) despite the perturbation ^V{t/T). 

Theorem. Consider Markovian dynamics satisfy- 
ing the above assumptions. In the hmit of large 
T, perfect adiabaticity is approached with an error 

that scales as 0(y^ Ts)' where A > is some 
fixed energy scale. In the adiabatic limit, the 
states inside 'H^i evolve under the unitary ?7^(s) — 
Texp (-z/; TVB{Pfery(g)P,r4l ® Qf{q)}dq). 

Proof. Let us divide the total time T into N steps, 
each of length 5t, T = NSt. We wiU take St = N/A 
(hence, T = iV^/A) such that when N — ^ oo, St is short 
on the time scale of change of the Lindbladian but long 
on the time scale for reaching the asymptotic regime 
of the instantaneous Lindbladian. The differentiability 
assumptions about C{s) and Vki{s) imply that we can 
write C{i^) = C{^) + 0(i), V^-^) = V{^) + O(^), 
t' e [Q,5t]. The evolution of the density matrix during a 
single time step can then be written 



9 r.A, 



p{t) ^ Pit + 6t) = Tefo' ^*'^(^)p(i) + 



(5) 



d,'e^^(*)(**-*') ( ^[y(l),e^~(^)*p(i)]) +0(^). 



Assume that the state at time t has the form 
Pit) = Ptiit) [Qfi^] + 0(1)] + O(^) 



(6) 



Then the first term on the right-hand side of Eq. ([5]) is 

Te/o" '^*'^(^)p(t) = p^^it) ® {gf{±) + 0(i)) + O(^), 
noiseless and Te-I^o -^(t^) = 



gStci^) ^ O(-l), so for large St the state on nf 
decays to ef (y) + 0{jf ) (see Appendix A for an exact 
relation to the decay rate) . For the second term, ignoring 
errors of order 0{^), we can use p{t) 



leaves p{t) invariant, so this term be 

rSt ,,, rl ±\(!\i-f'\lir/ t \ A/,\ ^ R 



But e^(^)*' 

comes ^ * di'e^(T)(**-t') [■(/(^), 
Using noiseless-subsystem properties of the Lindbla- 
dian 23, [i^l, in Appendix A we show that this term 
^J^'dt'e^<-j^i''-''^r,dVi^),pti{t) 
But C{s)Vki — 0, so the in- 



is equal to 



tegral yields ^Vm[V{^), pti{t) Qf{^)] = 

-zf [TrB{P,m^)^fer4l « efm^PkM ® efii) 

(the last inequality can be verified by a simple algebra). 

We therefore see that if the initial state is of the 
form ([6|), it will remain of this form for all times, up 



to an error 0(4 



^rp) resulting from the ac- 
cumulation of the errors O(j^) at every step. More- 
over, we see that the reduced density matrix on 
satisfies the difference equation p^i {t + St) — p^^ [t] = 

which in the limit N ^ oo yields the differential equation 



-AiiTs) = ~^[^rB{PCl"V{s)P^,''I^,®gf{s)},pUTs)] 
describing the effective evolution stated in the theorem. 

Note. Our theorem includes an adiabatic theorem for 
closed systems as a special case. However, the conver- 
gence rate stated in our theorem is weaker than the stan- 
dard one [the error is as opposed 0(^)] since 

our proof captures dissipative cases as well. (In Appendix 
A, we describe a natural energy scale A associated with 
the curve C{s), which can be regarded as a generalization 
of the minimum energy gap.) 

Decoherence- assisted computation in noiseless codes. — 
Computation in noiseless subsystems requires operations 
that keep the information inside the code 25|. However, 
the Hamiltonians that preserve the code in general may 
be rather complicated and may not be naturally avail- 
able in a particular experimental setup. Thus strategies 
for achieving encoded universality [26l | by other means 
are of particular interest [27| . An immediate implica- 
tion of the above theorem is that for the common case of 
time-homogenous Markovian noise with Lindbladian C 
[to play the role of C{t/T) in Eq. Q], any Hamiltonian 
perturbation ^V{t/T) acting during t € [0, T] would give 
rise to (possibly non-trivial) unitary evolutions inside the 
noiseless subsystems Tifj of £ within an arbitrary preci- 
sion for sufficiently large T. Thus given a set of available 
interactions {V^i} that can be turned on with variable 
strength, for a given subsystem Ti^i one can produce the 
set of effective interactions 



V, 



off 



Tr 



AB tA 



B{Pkl ^ti^kl ^kl 



of). 



(7) 



(Note that preparation of the states on Tif is not needed 
as they quickly decay to the fixed points.) Encoded uni- 
versality is achieved if the set { V^*^} spans the Lie algebra 
su{m) over "H^j. Remarkably this is possible even if the 
Hamiltonians {V^} commute (see example below). 

Such an approach was first proposed in Ref. [l^ for 
noiseless subspaces [AmiHf = 1) under certain noise 
models that can be interpreted as continuous Zcno mea- 
surements projecting onto the subspace. Equation ([55)1 
provides a generalization of this idea to noiseless subsys- 
tems (that may exist even when no noiseless subspaces 
exist) and arbitrary time-homogenous Markovian mod- 
els. As an example, in Appendix B we study a two-level 
noiseless subsystem of three spin-i particles under collec- 



tive decoherence The noiseless subsystem involves 



highly entangled states, and non-local interactions are in 
principle required to perform operations on the encoded 
qubit. However, we find that the decoherence process 
itself can be used to induce an effective universal set of 
gates on the code by acting with local Hamiltonians. 

Holonomic quantum computation via dissipation. — In 
the previous method, we assumed that the perturba- 
tion ^V{s) is applied by the experimenter. However, 
the conclusions are valid also if we assume that the de- 
scription is with respect to an instantaneous basis of 
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a time-dependent noiseless subsystem 'H^i{s) of C{s), 
where the perturbation now arises from the time depen- 
dence of the basis. As C{s) acts trivially on 'H"^i{s), the 
effective transformation in T-L^i^s) is not of dynamical 
origin. Indeed, in the adiabatic limit, an initial state 
p^^{Q) over 'H'^i{0) (XtHf (0) transforms via the superop- 
erator lim^^^o 'Pki{l)'Pki{l - 6s)...Vki{Ss)Vki{0) which is 
an intrinsically geometric quantity defined via the pro- 
jectors Vki{s). But the effective unitary on 'H^i{s) de- 
pends on the choice of basis for (s) and is not gauge 
invariant. However, if T-l'^i{s) is taken around a loop, 
H^iiO) = so that the final basis is the same as 

the initial one, the resultant transformation is a gauge- 
invariant quantity that generalizes the standard holon- 
omy associated with parallel transport of Hamiltonian 
eigenspaces We note that the idea of adiabatically 
"dragging" a subsystem (rather than a subspace) along 
suitable paths in order to perform geometric gates inside 
it has been proposed for the case of Hamiltonian dynam- 
ics as a powerful tool for robust computation [17|. How- 
ever, a subsystem cannot be dragged along an arbitrary 
path 'H'^(s) by a Hamiltonian since some paths necessar- 
ily give rise to correlations between T-L^{s) and 'H^{s). 
This problem does not exist here since the Lindbla- 
dian acting on 7i^{s) severs any such correlations. (For 
dissipation-driven holonomies in subspaces, see Ref. [28|.) 

The mathematical foundations of these geometric 
transformations will be studied elsewhere. Here we show 
that the method can be used for universal quantum com- 
putation. Consider a two-qubit system H — 
and a depolarizing Markovian channel acting locally on 
n^, ^ = C^p^ = - p^). Consider the unitary 

U{s) = e-*^('^'^^®'^?+^'^^''®^''), where V^?T¥ = 2n. 
The Hamiltonian in the exponent can be easily diago- 
nalized and one sees that U{1) = U{0) = I. Hence, 
if we change the Lindbladian via the unitary t/(s), we 
will take the noiseless subsystem V.^ around a loop 
■H^(s) (g) H^{s) = U{s)H^ (g) with a single- valued 
basis. According to Eq. (|36|) (here = ■^), the subsys- 
tem will experience the effective Hamiltonian bcr^ which 
gives rise to the transformation 6"*^*^= at the closing of 
the loop. Similarly, by exchanging az and cTx, we can 
generate the unitary e~'^'^^ . To perform an entangling 
gate between two qubits, A and A' , we can start with the 
same Lindbladian acting on B and rotate it via the uni- 
tary U{s) = e-*«("'^^'®'"^®'"?+'"^^'«"^^'»-f''), which gives 
rise to e~'^'^=" '^'^'^ . This set of gates is universal. 

Conclusion. — We introduced a theory of adiabatic 
Markovian dynamics that relates the notion of adiabatic- 
ity to the theory of noiseless subsystems. We proved an 
adiabatic theorem for such dynamics and proposed two 
novel methods of quantum information processing based 
on it — decoherence-assisted computation in noiseless sub- 
systems and dissipation-driven holonomic computation — 



that add to the developing picture of d issip ation as a 
powerful quantum computation primitive [29| . A natural 
problem for future research would be to find exact bounds 
on the adiabatic error in Markovian dynamics similar to 
those obtained for closed systems, e.g., in Ref. jsoj . 
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APPENDIX A: DETAILED PROOF OF THE 
ADIABATIC THEOREM 

Preliminaries 

Before we go in detail through the main steps of the 
proof, it is convenient to introduce an energy scale A 
associated with the curve C{s). This quantity can be 
regarded as a generalization of the minimal spectral gap 
of the Hamiltonian from the case of closed systems, and 
is a suitable choice in view of certain later calculations. 
Although it is not the purpose here, this energy scale 
could be useful for deriving exact bounds on the error 
and not just its scaling with T. 



As shown in Ref. 23 1, the subsystem is noiseless 



under the evolution driven by C{s), if and only if for every 
s the Hamiltonian and the Lindblad operators have the 
block forms 



His) 



/^®iff(.s) i/2(s)' 

i3,(s) 



(8) 
(9) 



where the upper- left block corresponds to H^i <E) H? , and 



(10) 



Then it is not difficult to verify (see also Ref. [2^) that 
C{s) preserves the subspace B2 of operators with vanish- 



ing lower right block, 

62 = {r e B{n) I r = 



T2 



T3 



(11) 



where B{7i) denotes the space of linear operators over 
v.. It further preserves the subspace of operators with 
vanishing lower right and offdiagonal blocks, 

61 = {r e B{n) I r 











where it acts as 



C{s) 





0' 






§£^(s)ri^^ 0" 
















(12) 



(13) 



where C^{s) is a local Lindbladian with Hamiltonian 
Hi{s) and Lindblad operators iy(s). Note that C^{s) 
has a non-degenerate eigenvalue with a corresponding 
right eigenoperator gf{s), and all its other eigenvalues 
have negative real parts, since gf (s) is an attractive fixed 
point. By continuity, the magnitudes of the real parts of 
the non-zero eigenvalues of C^{s) have a minimum value 
in the interval s e [0, 1]. Denote that value by Ai > 0. 

We will also need another quantity, A2, which is the 
minimum of the magnitudes of the non-zero eigenvalues 
of C{s)V2 in the interval s G [0, 1], where 7^2 is the pro- 
jector on B2 (this minimum also exist by the assumption 
of continuity in a closed interval). Then we can define 



A = min(Ai, A2), 



(14) 



which we will serve as a natural energy scale in our anal- 
ysis. If is a noiseless subspace {dirnHf = 1), Ai 
does not exist and A = A2. 

We note that in the case of closed systems, where 
is an eigenspace of the Hamiltonian, A is exactly the 
minimum gap that separates this eigenspace from the rest 
of the spectrum. This is because the non-zero eigenvalues 
of C{s)V2 are ±i{Ek{s) — Em{s)), m ^ k, where i?„ are 
the energies of the eigenspaces Hni- 

Remark. The existence of A > follows naturally 
from continuity and the assumption that dimH^; and 
dim are fixed during the closed interval s S [0,1]. 
The condition dimH^; — const(s) can be thought of as 
an analogue of the closed-system requirement that the 
eigenspace does not cross other energy levels. However, it 
may be possible to relax the condition dimHf = const (s) 
as long as we require A > for the open subintervals of 
s G [0, 1] during which dimHf = const(s). 

Before we proceed, we will need another observation. 
By the definition of "H^; ^Hf, the only right eigenopera- 
tors of >C(s) with eigenvalue inside the subspace B2 are 



those with T2 = and t 



AB 



^ki ® Qii^)- Denote the 



subspace of these operators by Bq{s) [Bq{s) C Si C B2Y 
Bo{s) = {r e B{H) I T 







(15) 
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Let ■P'(s) be the projector on the subspace B2 Q Bo{s). 
Then the superoperator £'{s) = C{s)V'{s) over B2 
Bq{s) has only non-zero eigenvalues and is therefore in- 
vertible. Its inverse, £'~^(y), has magnitude which is 
bounded by the inverse of A, 



c'-^ II < 



A' 



(16) 



The boundedness of this operator will be used at a certain 
stage of the proof, and the energy scale A was chosen as 
described above since it provides the bound ([TS]). 

We are now ready to go through the steps of the proof. 



Main proof 

Let us divide the total time T into N time steps, each 
of length 5t, T = N5t. We will take 5t = iV/A (hence, 
T = N^/A) such that when N ^ 00, St is short on 
the time scale of change of the Lidbladian but long on 
the time scale for reaching the asymptotic regime of the 
instantaneous Lindbladian. The evolution of the density 
matrix of the system during one time step can be written 

Pit) ^ Pit + St) = Tefo' + 



St 



dt"Tef'^^^"'^i^^ [^[Vii±f),. 
xre/o"*'^l^)p(i) + 0(-L), 



(17) 



where [V, -Jp — [V, p]. Our assumption that £(s) is 
differentiable and Vkiis) is twice-differentiable implies 
that £(s) and Vis) are differentiable. Hence we have 
Cit^) = £(±) + 0(i), Vii^) = Vi^) + 0(i), 
t' e [0, (5i]. We can therefore simplify Eq. (|T7)) : 



St 



pit) ^ pit + St) = Tefo' '^*'^(^V(i)+ (18) 
dt'e^m^*-^') (^[y(|),.]) e^~(*)*>(0 + O(^). 



Let us now assume that the state at time t has the 
form 



t . 1 



1 



Pit) = pUt) ® [sfi j;) + Oi-)] + O(^). (19) 
Then the first term on the right-hand side of Eq. ([T5|) is 
TeC ^ pA^t) ^ ,B^t) + (20) 

for some r^it), because "H^; is noiseless under £(s). We 
will now show that for sufficiently large St, 



rFit) = gfi^) + Oi^). 



(21) 



First of all, we have Te-io 



dt'C{i4j!-) 



eStci^)+Oi^ 



Also, under the action of Ci j^) (for fixed t) any state p 

gfi^), where p^ 



N I- 

AB 



on nil 



V.f decays towards p 



A 



TtbP^^ , with a rate at least A. Since St — N/A, we have 
e'5*^(T)p(0 =p^,(0«)[(?f (t) + 0(^)]. (In fact, the error 
Oij^) in the latter expression is an overestimate, but for 
our purposes this precision suffices.) Therefore, for the 
first term on the right-hand side of Eq. (fT8|) we obtain 

Telt ^^'^C-^)pit) = piit) ® [gfi^) + Oii)] + O(^). 

Next, consider the second term on the right-hand side 
of Eq. ((T5| . Ignoring terms of order 0(-^), for this term 
we can take pit) = p^iit) gfi^). The superoperator 
e-C(T)* leaves pit) invariant, so the expression becomes 



T 



St 



' dt'e^^^^^''~'')[Vi'-),pUt)®gfi'-)]. (22) 



We are now going to show that this expression equals 

(23) 

Indeed, let us add and subtract 'Pki\V i-j:), piiit)(^ gf i^)] 
from the operator \V i-^) , p^iit) ® gfi^)] in expression 
(I22t. We obtain 



St , , 

^ dt'e^ms'-'')r,,lVi^),pA^t)^gfil)] 



St 



^ / dt'e^<-r^'-'*-''^Wit), 



(24) 



where 



T 



Wit)^[Vi^),pt,it)^gfii)] 
-'PM[Vi^),ptiit)®gfi^)]. 



Note that [Vi j^),piiit)<g)gfij,)] belongs 10^62, and there- 
fore Wit) e B2Q Boi^) since Wit) decays to zero under 

the action of Ci-^), i.e., it has no component in Boij^). 
We can therefore formally solve 



St 



dt'e^(T)(«*-*')M/(t) 



.l,C'-Hj;){l-e^^^'^'')wit), (25) 



where £' ^i^) is the pseudo- inverse of £(|^) over B2 Q 
So(t)- According to Eq. (HH), || ||< 3, hence the 

magnitude of the term in Eq. (^5)) is Oi^^) = Oi^). 
The only non-trivial contribution to the expression (j24p 
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then comes from 

^^'*dt'e^~(^)(^-*-*')7',,[F(|),p^,(t)®,f(|)] (26) 

= -,^[TrB{F.r^(|)^.r4^ ® 9fi^)}, ptM ® ^^f (|) 

where in the first equality we used that C{s)Vki = 0, and 
the second equahty can be verified by a simple algebra. 

We now see that if we start with p(0) — p^iiO) (8) 
[pf (0) + 0{jj)], the state will remain of this form for all 

times, up to an error of order 0{jj) — ^) that re- 
sults from the accumulation of the errors 0{^) at every 
step (there are a total of TV steps). Moreover, the reduced 
density matrix in satisfies the difference equation 

Pi{t + 5t)-pti{t)= (27) 

which in the limit — > oo yields the differential equation 

^^pUTs) = -z[TrB{P,rV(.)F,r4l ® Qf{s)},pUTs)] 

(28) 

describing the effective evolution stated in the theorem. 
This completes the proof. 

APPENDIX B: EXAMPLE OF 
DECOHERENCE-ASSISTED COMPUTATION IN 
NOISELESS SUBSYSTEMS 



M = ®^l^i/2 ^nj ® M{dj), where J is the total angular 
momentum and M{dj) are dj x dj complex matrix alge- 
bras with multiplicity nj. In particular, dj = 2J-I- 1 and 
rij/2 2, 71,3/2 = 1. The Hilbert space correspondingly 
decomposes as 

71^%^®%^®%^, (30) 

where 'H'^ is a noiseless qubit subsystem arising from the 
two- fold multiplicity of A^(2), T-L^ is the noisy cofactor 
supporting (2), and 'H'" is a noisy subspace supporting 
M.{A). The subsystems 'H^^'H^ can be described in the 
basis 



|0)-4|0)^ 


= -L(ioii)-iioi)), 




(31) 


|0)^|1)^ 


- i=(|010)-|100)). 




(32) 


|1)^|0)^ 


= -^(2|1I0)-|101)- 


1011)), 


(33) 




= -L(_2|001) + 1010) 


+ I100)). 


(34) 



One can verify that under the evolution (P^ . there is a 
unique fixed point on "H^, 

q"" - — ^I0)(0|^ + — ^|1)(1|^. (35) 
7-1-7^ 7-1-7^ 

Similarly, there is a unique fixed point on the subspace 
n^. Using that 



To illustrate the idea of decoherence-assisted quantum 
computation in noiseless subsystems, we consider as an 
example a two-level noiseless subsystem of three spin-i 
particles under collective decoherence 14| . 

Under the evolution [sij 

^ = -iuj[J,_,p\+-i~{J-pJ+ - ^J+J-p- ^pJ+J-) 

+j+{J+pJ- - \j-J+P - \pJ-J+). (29) 

where Jz — h'^l J± — ^ ■ a\. are collective spin 
operators and 7^, 7^ > 0, there are no non-trivial noise- 
less subspaces but there is a noiseless subsystem. The (f- 
closed) operator algebra generated by J a is isomorphic to 



v: 



off 



AB 
ij 



(36) 



we obtain that the local Hamiltonian a\ gives rise to 
the effective Hamiltonian — 7-, ~^ — ttCf^, where — 

2V3(7~-|-7+) ^ ^ 

|0)(1|^ + |1)(0|'^ is the encoded Pauh operator on "H^. 



Similarly, a\ gives rise to 



2(7"-7+) 



cr^ . These two Hamil- 



3(7-+7+)^ 

tonians generate SU{2) on H^. 

For universal computation one needs the ability to en- 
tangle multiple noiseless qubits, e.g., by bringing different 
blocks together [iH and manipulating the logical infor- 
mation inside the resulting larger noiseless subsystems. 
This problem can be treated via the same approach. 



